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We calculate various properties of pseudoscalar mesons in partially 
quenched QCD using chiral perturbation theory through next-to-leading or- 
der. Our results can be used to extrapolate to QCD from partially quenched 
simulations, as long as the latter use three light dynamical quarks. In other 
words, one can use unphysical simulations to extract physical quantities — in 
this case the quark masses, meson decay constants, and the Gasser-Leutwyler 
parameters L4 — Lg. Our proposal for determining Lj makes explicit use of an 
unphysical (yet measurable) effect of partially quenched theories, namely the 
double-pole that appears in certain two-point correlation functions. Most of 
our calculations are done for sea quarks having up to three different masses, 
except for our result for L7, which is derived for degenerate sea quarks. 
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I. INTRODUCTION 



A major obstacle to direct simulations of lattice QCD is the difficulty in simulating 
with light dynamical quarks. In particular, the up and down quarks must be reached by 
a chiral extrapolation. In present simulations one must do this extrapolation from roughly 
ms/2, where is the physical strange quark mass. This is far from the light quark masses 
(m = (m„ + md)/2 ^ ms/25). 

The aim of this paper is to provide formulae which can aid in this extrapolation. To 
do this we use chiral perturbation theory (ChPT) at next-to-leading order (NLO). The 
parameters of the chiral Lagrangian that enter at this order are the Gasser-Leutwyler (GL) 
coefficients, Li—Liq. An alternative way to view the extrapolation to QCD is that, by ffiting 
numerical results in a region where quark masses are considerably larger than the physical 
light quarks, but small enough that NLO chiral perturbation theory is sufficiently accurate, 
one determines the relevant Lj. These are physical parameters of QCD, governing many 
different physical properties (e.g. pion masses and scattering amplitudes). With the Lj in 
hand, one can then extrapolate to QCD, and, in particular, determine the physical light 
quark masses. For example, as has been stressed in Refs. determining the combination 
2^8 — L5 with only moderate accuracy might allow one to rule out the interesting possibility 
that rriu = 0. The accuracy of extrapolation depends, of course, on the reliability of NLO 
chiral perturbation theory. This can be studied by seeing how well the numerical data fit 
the expected forms, including the curvature predicted by chiral logarithms. 
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FIG. 1. Schematic representation of the space of PQ theories for "light" quarks (defined as 
lighter than the physical strange quark mass). The approximate range of present simulations is 
shown. The shape of this region is determined by the fact that critical slowing down is less severe 
for valence quarks than for sea quarks. 

An observation of practical importance is that one can make use of partially quenched 
(PQ) simulations to aid in the extrapolation to QCD In partially quenched simulations, 
one changes the mass of the valence quarks (typically reducing them), while holding the 
dynamical (or "sea" ) quark masses fixed. The situation is illustrated schematically in Fig. |l|. 
This leads one into a space of unphysical theories, from which one might expect to obtain only 
qualitative information about QCD. It turns out, however, that, if all quark masses (valence 
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and sea) are small enough, one can use PQ theories to obtain quantitative information about 
unquenched theories. Since it is computationally less demanding to reduce valence quark 
masses, PQ simulations are often used to obtain approximate information on QCD. Our 
point here is that they can be used to obtain exact information about QCD. 

This observation follows from the structure of chiral perturbation theory (ChPT) gen- 
eralized to partially quenched theories — PQChPT 0. The key point is that there is a 
subspace of quark masses (corresponding to the diagonal line in Fig. |I|) where PQChPT 
is completely equivalent to chiral perturbation theory for unquenched, QCD-like theories. 
Since the quark mass dependence in PQChPT is explicit, as in ChPT, it follows that the 
parameters of the PQ chiral Lagrangian (with 3 light sea quarks) are the same as those 
of QCD. These parameters do, however, depend on the number of sea quarks, N . This 
means that PQ simulations with = 3 light sea quarks, whatever their precise masses, give 
information about the parameters of the chiral Lagrangian of QCD. On the other hand PQ 
simulations with = 2 or iV = 4, the values used in most simulations to date, do not give 
direct information about QCD, even after extrapolation. 

These comments motivate the calculation of the NLO chiral corrections to physically 
interesting quantities in PQ theories. Some results already exist in the literature: those for 
charged pion masses and decay constants with degenerate sea quarks heavy-light meson 
properties 0, vector and tensor meson properties 0, baryons masses at large A^^^ 0, and 
electroweak amplitudes 0. We provide here two new results. First, charged pion masses and 
decay constants are considered for non-degenerate sea quarks (having up to three different 
masses). This completes the calculations of simple pion properties for any theory that is 
likely to be simulated. It allows one to determine L4_6 and Lg. Note that, although non- 
degenerate sea quarks are not necessary in order to extract these GL coefficients, as has been 
stressed in Refs. there is no drawback to using them. Indeed, someone might prefer 

to extrapolate using more "QCD-like" simulations with two degenerate "light" dynamical 
quarks and one dynamical quark with mass fixed close to the physical strange quark mass. 
Our formulae apply for such a theory. 

Our second new result concerns the GL coefficient Lj. In QCD, this contributes only to 
neutral meson masses, and does so proportional to (m^ — m)^. To determine L7 using meson 
masses from unquenched simulations thus requires non-degenerate quarks. At first sight, 
PQ simulations do not improve the situation: neutral meson masses are still independent 
of L-j when the sea quarks are degenerate. We find, however, that one can determine L7 
from the coefficient of the double-pole in the propagators of neutral valence mesons, even 
for degenerate sea quarks. This is a nice example of the utility of PQ theories. Although the 
double-pole is itself an indicator that these theories are unphysical, its effects can nevertheless 
be measured in lattice simulations, and its inferred coefficient turns out to be related to a 
physical quantity. 

Throughout our calculations, we treat the rj' as a heavy particle, and integrate it out, 
following Ref. @]. This greatly simplifies the resulting expressions, since it removes the 
dependence on additional rj' coupling constants. It raises, however, two important issues. 
First, is the 77' heavy enough that removing it is appropriate? The answer depends on the 
number of light sea quarks, A^, and the number of colors, N^- For the physical values of 
these parameters, A^ = A'c = 3, we know from experiment that Mr^i ?s 1 GeV. Since this 
is the scale at which chiral perturbation theory breaks down, it is appropriate to integrate 
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out the Vj for these theories. We stress, however, that our formulae are only valid when 
the dynamical quarks are light enough that all pseudo-Goldstone bosons, including those 
composed only of sea quarks, satisfy MpQ^ <C M^,. For further discussion of this point, and 
of the limitations of the approach taken here, see Ref. f^. 

The second issue is more technical. How does one integrate out the rj' in PQ theories? 
In this paper we follow Ref. |Q and do this by hand, working only at tree level. This is 
unsatisfactory, since, as we know from QCD, rj' loops to all orders give contributions of the 
same order in chiral perturbation theory. In other words, the rj' must be integrated out 
non-perturbatively. We return to this issue in a companion paper, where we demonstrate 
that the approach adopted here is equivalent to integrating out the r]' non-perturbatively in 
the PQ theory ^. 

This paper is organized as follows. In the following section we recall the formalism of 
PQ chiral perturbation theory, and explain our calculational procedure. After presenting a 
simple form for the neutral meson propagator in Sec. [TTTI , we give our results for charged 
pion properties in Sec. We analyze these results in Sec. 0, paying particular attention 
to the convergence of the chiral expansion and the importance of non-analytic terms. In 
Sec. 0we explain how to extract L7 using PQ theories. We end with some conclusions. Two 
appendices deal with technical issues in the calculation of the neutral meson propagator. 

Some parts of this work have been reported previously in Ref. M . 



II. THEORETICAL FRAMEWORK 

We consider partially quenched theories with the following quark complement: A^i + 
N2 + N^i = N sea quarks, Ni each of mass mf, two valence quarks with masses and m^; 
and two corresponding ghost quarks with masses = mA and = tub- The ghosts are 



needed to cancel the determinant arising from the valence quark functional integral |jTO[ . In 
the chiral limit, this theory has an SU{N + 2\2)l ® SU{N + 2|2)/j symmetry group If 
Ni = N2 = = 1, and mi = m2 = and = iris, then the sea-quark sector is 
QCD. Generalizing to arbitrary numbers Ni covers most other theories that are likely to be 
simulated in an effort to shed light on QCD. 

An important property of PQ theories, which follows trivially from their definition, is 
that the sea-quark sector decouples from the valence sector. To be precise, all correlation 
functions composed of only sea-quark fields are identical to those in the unquenched sea- 
quark theory. There is no "back-reaction" from the valence sector. The same result must also 
hold for the low energy chiral Lagrangian describing the PQ theory: correlators of pseudo- 
Goldstone mesons composed of sea quarks should be the same as in the chiral Lagrangian 
describing the unquenched theory. This was shown to be true in Ref. 0. In practice, 
however, one might view all correlation functions that are calculated as being those of 
valence quarks, and so it is more useful to reformulate this property as follows. When 
each of the valence quarks is assigned a mass that is equal to one of the sea quark masses, 
sea and valence quarks become indistinguishable, and the cancelations of the "doubled" 
quark species against their ghost counterparts trivially render the theory the same as an 
unquenched theory containing only the "original" sea quarks. That this is so was also 
shown in Ref. P], and it has important consequences in the following. 
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At low energies, the partially quenched chiral effective theory is expressed in terms of 
the fields HQ 



i'^ab) = ^ (nab) 



$ = 

$0 = str$/ViV, 
and the quantity 

X = 2fim = 2/i diag(mA, mB, mi, 



a,be {A,B,1,2,...,N,A,B} 



mi, 1712, 



m2,m^. 



m3,mA,mB 



(2.1) 
(2.2) 

(2.3) 
(2.4) 



N2 



N3 



where m is the quark mass matrix. In the following we use the notation xa = '^fJ'^A, 
Xi = 2/imi, etc. The constants / and fi are unknown parameters. The fields iiab describe 
the pseudo-Goldstone particles of the theory — we refer to them generically as mesons even 
though some are fermionic. Because of the anomaly, arbitrary functions of the field $o (the 
super-?]') can appear in the Lagrangian. 

The partially quenched chiral effective Lagrangian is expanded in powers of ~ 
M^/A^ ~ / h?, where M is a typical pseudoscalar meson mass, p the momentum, and 
A ~ 1 GeV is the scale beyond which the theory breaks down. 

The parts of the Euclidean Lagrangian contributing to meson masses and decay constants 
at one-loop order are: 



str (d^W^T)) - ^str (^S^ + Sx) + ad^^^d^^^ + ml^l 
U str (a^Sa^St) str (xSt + Ex) + str [S^SS^S^ (^S^ + 



NLO,2 



Str {xT.^ + Sx)] - str (xS^xS^" + Sx^x 
str (xSt - Sx)l' + vi^l^ii (S^Sa^St 



-L' 



+ V2%sii[xT.' + T.X 



(2.5) 



(2.6) 



(2.7) 



The coefficients a, uiq, Lj and f j are further unknown parameters of the low energy theory.^ 
The Li depend, in general, on the renormalization scale. The NLO Lagrangian is broken 
into two parts because flavor off-diagonal mesons receive contributions only from £nlo,i- 

At this point we can make clear the relationship between the partially quenched chiral 
Lagrangian and that describing low energy QCD. The latter is obtained by setting = 3, 
and "unquenching" — i.e. assigning and itlb values from {mi,m2,m3}. It follows that 
the unknown coefficients in C are, for A^ = 3, identical to those in the QCD chiral Lagrangian. 
This shows that these constants also govern the chiral behavior of PQ extensions of QCD. 



^The normalization of is different from that used in Ref. [P, although it agrees when A = 3 
as in QCD. 



We revert to the notation a of Ref. [ 11 1 , rather than the used in Ref. @ . 
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In QCD, one can take a further step and "integrate out" the i]' . This is appropriate since 
it is not a pseudo-Goldstone boson, having M^, pa itIq + 0{m) 1 GeV^. Technically, the 
matching between theories with and without the r]' is non-perturbative, since loops involving 
the Vj are not suppressed by powers of M ^ or p^. Thus in the standard approach one simply 
writes down the Lagrangian without the r/', and it has the same form as Eqs. ( p.5| )-( p77D , 
except that $ is traceless. It follows that a, mo and the Vi are irrelevant, and the only 
NLO coefficients are the Lj. It is in fact in this theory that the — the Gasser-Leutwyler 
coefficients — are conventionally defined. 

In previous work on PQQCD, the step of integrating out the r/' has been done by hand, 
i.e. at the level of individual diagrams rather than the Lagrangian [Q. We summarize the 
procedure here — details will become apparent in the following section. 

• Loop diagrams involving the v( are dropped, since these lead to shifts in the parameters 
Li which are automatically included if we use the Lj from the QCD Lagrangian without 
the rj . 

• Couplings special to the rl , such as the Vi in Eq. ( p.7| ), are treated as small, of O(e^), 
and thus appear only at tree level. The justification for this treatment is that these 
couplings are suppressed by powers of l/Ai'c, in this case l/N"^- 

• On the other hand, the parameter is treated non-perturbatively since it is known 
to be ~ A^, despite the fact that it is proportional to 1/Nc. In particular, we treat 
M'^/ml as O(e^) (with M, as above, a typical meson mass). For convenience, we also 
treat a non-perturbatively. 

While this procedure may be accurate enough for phenomenological purposes, it is theoreti- 
cally unsatisfactory because the rj' should be integrated out non-perturbatively. As noted in 
the introduction, we will address this concern in a separate paper 0. In particular, we will 
argue that the procedure adopted here in fact leads to results that are equivalent to those 
obtained from integrating out the rj' non-perturbatively. 



We close this section by deriving a result needed in Sec. |V|. We claimed above that 
discarding rj' loops allows us to write our results in terms of the standard GL coefficients, 
i.e. those in an effective Lagrangian without the rj'. This is not quite correct. Tree diagrams 
involving an intermediate t]', which we keep, lead to a shift in L-j proportional to M^/rriQ. 
Thus we are not using the conventional Lj, but rather that defined in the effective theory 
containing the t]', which we denote L'^ [as anticipated in Eq. ( p.7|) ]. In order to express 
our results in terms of conventional parameters, we need to relate Lj to L7 within the 
approximations of our procedure. 

To determine this relation we need consider only the sea-quark sector, and thus work 
with the conventional (unquenched) chiral Lagrangian. The t]' field in this theory is the 
restriction to the sea sector of the "super-?]'" field $o- To make the t]' dependence explicit, 
we decompose E into pseudo-Goldstone and rj' parts: 

S = t/exp(-^); UeSU{N), (2.8) 

and substitute into the chiral Lagrangian. The result is the original form with —>■ U and 
$0 ~^ 0, i.e. the standard QCD chiral Lagrangian, plus the following terms 
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= Cl<l>0 + C2^l + C3(a$o)' + 0{^l) (2.9) 

Cl = iMxU^ - xU)/{2Vn) + O(x') (2.10) 

C2 = m^ + 0(x) (2.11) 

C3 = l + a + 0(x) (2.12) 

Only the leading order terms in the chiral expansion of the coefficients are shown, since 
higher order terms give contributions to the conventional chiral Lagrangian of orders and 
higher, too high to effect the matching of the GL coefficients which appear at order e^. For 
the same reason, the (d^o)^ term can be dropped. Keeping the rj' in tree graphs amounts to 
doing the functional integral over $o keeping only linear and quadratic terms. The result is 
a contribution to the conventional chiral Lagrangian of order with the same form as the 
L7 term: 

MxU-xU^)?{l+Oix))- (2.13) 



16Nml 

Thus we find, within our approximations, the relation 



f2 

Lt = L' — ^. (2.14) 



III. CALCULATION 

Expanding Llo to quadratic order, we obtain the meson propagators. For "charged" 
mesons (i.e. flavor off-diagonal states) these are 

G'iM = / rf^xe-^^-(7r,,(x)7r,.(0)) = ^^^^2^x^)1'^ ^ ' ^^'^^ 



where the signature vector is 



1 a e {A, 5, 1,2, 3,..., AT} 
-1 a G 



(3.2) 



The propagators for "neutral" (flavor-diagonal) mesons include the contributions of the 
super- 77' interactions. A general expression has been given in Ref. [Q], but we find it conve- 
nient to use an alternative form. The propagator is a matrix acting on the space of neutral 
meson fields, i^aa-, a = 1, + 4. In app. |^ we show that 



(3.3) 



^ ea6ab _ {ml + ap^)/N {p^ + xi) {p^ + X2) (p^ + X3) 

+ Xa {p^ + Xa) (p2 + Xb) (1 + a) (p2 + Ml) (p2 + M^) (p2 + M^, ) ■ 



Here the ttq, t] and 77' are neutral mesons in the sea-quark sector. For = 3 they are 
the usual neutral mesons of QCD; for > 3 they are the appropriate generalizations, as 
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explained in the appendix. Their masses are functions of the sea-quark masses, and of the 
Ni, as given exphcitly in Eqs. ( |A19| )- (|A21|) . 

The neutral propagator shows explicitly the unphysical nature of the PQ theory. For 
example, if a = b = A, the second term has a double-pole at = — Xa- These double-poles 
are absent, however, in the physical, sea-quark, sector. This was shown in Ref. 0, but 
is particularly transparent with our result. For example if a = 6 = 1 (or equivalently if 
a = b = A and xa = Xi)? then the (p^ + Xi) i^i the numerator reduces the double-pole to a 
single-pole. 

For calculations, it is preferable to rewrite the propagator as a sum of (single or double) 
poles. For simplicity, we discuss the case when a ^ b and Xa 7^ Xfe? fo^' then has only 
single poles. For uniformity of notation, we introduce the definitions 

X. = Ml^, X, = M^, X,'=M^,. (3.4) 

in terms of which 

GS = ~ E ^ i-^"), (3.5) 

{ml - axx) Y{i=ifi{Xt - Xx) 
[I + a)Y{y^^{Xy - Xx) 

where x and y run over a, 6, vr, 77, rj' . The degenerate limit Xa = Xb 01 a = b, which, in 
general, has double poles is straightforward to obtain. 

At this point we are in a position to integrate out the rj' by hand, bearing in mind that 
this propagator appears in loops. First, as discussed above, we drop the rj' pole. Second, 
we expand the residues of the remaining poles in powers of Xx/^o, x 7^ rj' , and drop all 
but the leading term. This is justified since we use the neutral propagator in one- loop 
diagrams, which already give NLO contributions. The only exception is in our discussion of 
Lj in Sec. |VI|, where the neutral propagator appears at tree-level. With these changes, the 
residues become 

nj=i,3(x« Xx) /„ _,s 

Rx = — 7 T- {i-l) 

[{v^a.b.^,ri[Xy - Xx) 

Note that both mg and a disappear from the neutral propagator. 



IV. NLO RESULTS 

In this section we calculate the properties of a meson composed of two valence quarks to 
NLO. We use dimensional regularization, and subtract the poles following the conventions 



of Ref. 12 



Its mass, Mab, is obtained from the diagrams of Fig. ^. We find 

Ml^ = ^4^(l + Ce + Cp) (4.1) 



2 

''tree y2 K^^i) -^ij A i y2 



^ - ^ {2L, - L4) X + « {2Ls - L,) ixA + Xb) (4.2) 
1 



S\oop = iQp^2p^ {RaXa log Xa + RbXb log Xb + RnXn log Xtt + Rr,Xv l^S Xv} (4-3) 
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Here x is the average sea-quark mass, 



X=^JT.N^X^. (4.4) 
i=l,3 



and the residues are defined in Eq. ( p.7D . For concreteness we quote two exphcit examples 

{Xa - Xi) {Xa - X2) {xa - X3) 



Ra 

Rtt 



{Xa - Xb) {xa - X-k) {xa - X-n) 
(Xtt - Xi) (Xtt - X2) (Xtt - X3) 
(Xtt - Xa) (Xtt - Xb) (Xtt - Xr,) ' 



Rb is obtained from R^ by interchanging A and i?, and i?^ is obtained from R^^ by inter- 
changing TT and 7]. 



AB AB AB__0__AB 

NLO 

FIG. 2. Diagrams contributing to Mab- The letters next to the hnes denote the flavor indices 
of the propagating mesons. "VV" stands for Valence- Valence meson with V = A, B. "LO" and 
"NLO" describe the order of the vertex that makes the diagram contribute at 1 loop. 

The renormaUzation scale is implicit in the logarithms and the L,. Using the result 

E X.R.= {XA + XB-x) + 0{xVml) (4.5) 

X=A,B,1T,T] 

we see that a change in renormalization scale can be absorbed by shifting the L^. We 
have checked (for N = 3) that the scale dependence of the Li in QCD does render M^^ 
independent of the renormalization scale. 




VV, VS 



AB AB AB _ AB AB V / AB 
BWWWA 




LO, NLO NLO LO LO 



AB 



VV, VS 
AB 




LO LO 

FIG. 3. diagrams contributing to Jab- The wavy line represents the insertion of the axial current 
operator jt!^ABiP)- "^S" stands for Valence-Sea mesons with V = A,B and 5 = 1,2, 3. 
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To determine the meson decay constant, /ab, the axial current j^^^ is calculated to 
NLO, and then we evaluate the matrix element 



{0\j^AB\^AB{p))=tV^fABp' 



(4.6) 



using the diagrams of Fig. ^. The result is 



Jab — / (l + ^Lcc + ^YS loop + '^VV loop 
A AT O 



C>tree = -J^X^A. + J^{XA + Xb) 
jrf Ar ^ + fXA + Xi\ , 



•^VS loop 



=1,3 
1 1 



B) 



(4.7) 
(4.8) 

(4.9) 



VV loop 4^ ^g^2 J2 

[Xa-Xb) l 

, / D / N riog(x^/XA) log(x^/x. 

+ X-kRAXb-Xa) 

\ V Xa-X-k 



(4.10) 




where 



Y\i=i,Axi - xa) 

(Xtt - Xa){Xv - Xa) 



(4.11) 



and Db obtained by [A ^ B), are the coefficients of the double poles in the neutral prop- 
agators. It is straightforward to see that the scale dependence can be absorbed by shifts in 
L4 and L5, and we have checked that these shifts are consistent with standard results for 
AT = 3. 

As noted in the introduction, our formulae can be used to extract the GL coefficients by 
fitting to results from simulations. PQ simulations allow one to vary the valence and sea- 
quark masses independently, and thus to separately determine La^_q and Lg. In fact, since 
the NLO analytic terms, Eqs. ([4.2|) and ( [4.8|) , depend on the quark masses only through the 
combinations xa + Xb and x, one need only consider degenerate sea quarks.Q Indeed, to 
extract 2L^ — L5, the combination which determines whether m„ = 0, it is sufficient to use 
a single sea quark mass (as long as it is light enough that the formulae apply), and vary the 
valence quark masses. By contrast, with unquenched simulations, one would have to use 
non-degenerate sea quarks to separately determine all four L's. 

The most likely practical applications of these results are for simulations done with 2 
rather than 3 types of nondegenerate sea quarks. For QCD this would correspond to the 
limit of exact isospin symmetry. The results for this case can be obtained from those given 



^Note that the NLO analytic dependence of is not the most general quadratic term, sym- 

metric under A <-> i?, and vanishing when XA = XB = 0. Such a form would contain a term 
proportional to {xA — Xb)"^-, which is in fact forbidden by chiral symmetry. 
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above by carefully taking the limit X2 Xi- When two types of quark are degenerate, 
one of the neutral sea-quark eigenstates becomes an exact flavor non-singlet, and we choose 
this to be the pion. Thus we have Xn = Xi = X2- The remaining neutral state has mass 
Xri = Xi + X3 ~ X- III this limit, the residues simplify, e.g. 

P iXA - Xi) iXA - X3) 

[Xa - Xb> [Xa - Xr,) 
— ^ (4.13) 

^ (Xr? ~ Xl) {Xy ~ Xs) j^^-j 

iXv - Xa) iXv - Xb) 

with a similar cancelation in the Da,b- With these changes, the formulae given above still 
hold. 

We have checked our results in two ways. First, if we take all sea quarks to be degenerate, 
we obtain the results of Ref. 0.0 Second, we can consider the unquenched limit. By choosing 
Xa and xb to be equal to combinations of Xi-3 ^6 obtain the correct one-loop form for the 
masses of the vr"'", and . 

The residues R^ and are singular when pairs of the x's become degenerate, e.g. xa 
Xb and xa Xn- As expected, however, these singularities cancel in the full expressions 
for Mab and J'ab, which are analytic functions of the quark masses except in the massless 
limit. 

In Ref. 01 , it was emphasized that the one- loop corrections can diverge when the valence- 
quark masses are sent to zero at fixed sea-quark mass, leading to a breakdown of chiral 
perturbation theory. This discussion was based on the results for degenerate sea quarks, 
and we can now generalize it to non-degenerate sea-quarks. For the meson masses, the 
possibly divergent contribution is S^^^^, and we see that this only diverges if both niA and 
rriB vanish, in fixed ratio, but not if only one vanishes. For the decay constant the pattern 
is opposite: ^yy ^oop diverges if one of the valence masses vanishes, but not if both do in 
fixed ratio. This is the same pattern of divergences as for degenerate sea quarks; the sea 
quark masses only influence the coefficients of the divergent log xa and log xb terms. In the 
following section we discuss the practical implications of these divergences. 

It was also noted in Ref. that one can form combinations of the squared meson 
masses and decay constants from which the analytic correction terms (^tree) cancel. One 
thus predicts these combinations in terms of the quark masses and the leading order chiral 
coefficients fi and /, up to NNLO corrections. Studying these combinations in simulations 
allows one to test the applicability of NLO chiral perturbation theory. What we want to 
point out here is that exactly the same quantities can be used for non-degenerate sea-quarks: 
the Li still cancel. We do not, however, give the explicit expressions since they are lengthy 
and unilluminating. 

Finally, we discuss to what extent our formulae can be applied to lattice results obtained 
at non-zero lattice spacing. For definiteness, we first consider a calculation using Wilson 



^Except for the following typos in Ref. pointed out by Jochen Heitger, Rainer Sommer and 
Hartmut Wittig: in Eqs. (18), (19) and (20) 04 should be replaced by 04/2. 
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fermions, in which we work at a fixed bare couphng and vary the bare valence and sea quark 
masses. Meson properties are of course calculated in lattice units, i.e. we obtain M^^a and 

There are two effects of working at a non-zero lattice spacing. The first is that chiral 
symmetry, upon which our calculation is based, is broken explicitly. This symmetry breaking 



can, however, be incorporated into the chiral Lagrangian framework, as shown in Ref. [jT3|. 
The result is that all effects of the explicit symmetry breaking are of 0(a), except for the 
additive renormalization of the quark masses. Because of this, our formulae are valid, up to 
corrections of 0(a), as long as one uses so-called vector or axial Ward identity quark masses. 

Note that the corrections of 0(a) cannot be incorporated into our formulae by simply 
introducing a dependence into the parameters of the chiral Lagrangian parameters, /, /i and 
the Li. There are additional unknown constants which enter. 

The second effect is that the lattice spacing itself depends on the quark masses, at fixed 
bare couphng. This introduces an additional mass dependence into quantities expressed in 
lattice units. We note, however, that the mass dependence of a is a discretization effect of 
0{am) induced by the explicit chiral symmetry breaking |T^. Indeed, for non-perturbatively 
on-shell-improved Wilson fermions one can calculate this dependence, with the result 

a{m) = a{m = 0)[l-cam + ...], c = -^^^^ , (4.15) 
where m is the average sea-quark mass, bg{g) oc g"^ + 0{g^) is an improvement coefficient 



introduced in Ref. |]T4|, and Plat = —dg/d\na oc g^ + 0{g^) is the lattice /3-function. Note 



that this 0(a) effect can be absorbed by shifting the GL coefficients as follows: 



(4.16) 



Alternatively, one could adjust the bare coupling as the quark masses are varied so as to 
keep the lattice spacing fixed. 

In summary, our formulae are approximately valid for meson properties expressed in 
lattice units, with the errors being of 0(a). Some, but not all, of these discretization 
errors can be absorbed into a dependence of the parameters of the chiral Lagrangian. With 
staggered or overlap fermions the errors would instead be of O(a^). Since discretization 
errors can still be substantial at present lattice spacings, it may be better to extrapolate 
first to a = 0, and then fit to the predicted forms. 



V. BEHAVIOR OF THE CHIRAL EXPANSION 



In the framework of chiral perturbation theory (regardless of quenching) one assumes 
that for any quantity calculated to a given order in the chiral expansion, higher order terms 
are smaller. Once the unknown couplings and parameters of the theory are determined 
(e.g. by experiment, or by fitting to lattice data) the consistency of the computation can be 
checked numerically. This is the main purpose of the current section. 

We can make this check because, using experimental data, we have reasonable estimates 
for the actual values of some of the Gasser-Leutwyler coefficients. Thus we can use our 
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results to predict the meson masses and decay constants in PQ simulations. Of course, these 
predictions are approximate because we only know the Lj approximately — this is where the 
lattice results themselves come in — but they give us a reasonable idea of how the chiral 
expansion behaves. 

We consider QCD with exact isospin symmetry, i.e. A^^i = A^2 = ^3 = 1 and mi = 
^2 = rriu and 1713 = rris. The meson masses and decay constants then depend on the 
seven parameters /, Xu, Xs, -^^4-6 and Lg. We want to choose these parameters so that 
the four charged meson quantities M^, M^, fn and fx take their experimental values. In 
order to match the number of parameters and observables, we take as starting values the 
GL coefficients quoted in Ref. [|T^. These are based on experimental results {L^{M,^) 
2.3 10-^ L^{M^) ^ 1.2 10-3) and the large limit {Li{M^) ^ Lq{M^) ^ 0). We then 
take our four free parameters to be /, Xs and the scale, A^, at which Li take the values 
just quoted. In effect, this moves us through the space of Lj on a particular path, which, 
when Kl M^, is consistent with our knowledge about the Lj. We claim no fundamental 
basis for this path — we use it for simplicity. It allows us to determine the dimensionless 
quantities //Al, Xu/ P and Xs/ P by fitting the ratios M^/fx, Mx/fx and fK/f-K, and 
then to determine A^ by requiring, say, M^^ = 140 MeV. We find A^ = M^(l — .0021) (using 
Mjj = 547 MeV), so that the fitted Lj's are, in fact, close to the inputs. The other outputs 
are / = 85MeV, Xu/Xs = 0.044 and Xs = (673MeV)2. 

We stress that we are not claiming that we have found a unique set of parameters. There 
is a region in the space of the Li which can describe the experimental observables, and for 
which it turns out that the chiral expansion is under reasonable control. We have picked, 
somewhat arbitrarily, one point in this region. 

We can now explore the behavior of the chiral expansion as a function of the four quan- 
tities 

Xa Xb Xi Xa 
7 5 5 

Xs Xs Xs Xs 

We have chosen to normalize the various quark masses relative to the physical strange quark 
mass, so that a ratio of unity represents the outer limit of where one would expect chiral 
perturbation theory to be reliable. We consider two types of two-dimensional cross-section 
of the parameter space: {y, y, x, 1) and {y, 1, x, 1). In both cases y corresponds to a valence 
mass while x is proportional to a sea quark mass. We name these cross-sections the "vr-plane" 
and the "i^-plane" because the unquenched line y = x in the former describes a pion-like 
meson made up of two identical light quarks, whereas, in the latter, y = x corresponds to a 
kaon-like meson. 

We examine the relative size of the NLO contributions by plotting 

(^L)NLO/(^li^)LO = (Co + Cp) = (5.2) 

and 

(/^b)NLO/ (/^s)LO = '^toee + 4s loop + 4v loop = (5-3) 

[see Eqs. ([4.1| ), ([4.7|) ]. In each plane, these functions are plotted along rays emerging from 
the origin at angles 15°, 30°, 45°, 60° and 75° with respect to the x-axis, confined to the 
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unit square (Fig. ^). The 45° line corresponds to unquenched theories, the other hues to PQ 
theories. The plots are shown in Figs. We also show a contour plot of 5-^ in the K-plane 
in Fig. |. 
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0.2 0.4 0.6 0.8 1 

FIG. 4. The dashed lines are rays in the 7r-plane or if-plane along which and are plotted 
in Figs. HH. 



5^(j = jctanG,^) 




0.2 0.4 0.6 0.8 1 



x = — 



FIG. 5. b^ in the vr-plane, 1), is plotted along rays of angle Q with respect to the x-axis 

(see Fig. Q). The values of theta in degrees are indicated next to the corresponding curves. The 
point that corresponds to the "physical" pion, iXuiXuiXu^Xs) is labeled "QCD". 
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5^(y = xtan 0,x) 
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5;^(3; = xtane,x) 




FIG. 8. Plots of 5^ in the iC-plane. 




0.2 0.4 0.6 0.8 1 



FIG. 9. Contours of (5^. The values of 5^ range from < .275 in the black region on the right 
up to values > .5 in the white region. There is a difference of .025 between the values of 5^ on 
adjacent contours 

The first conclusion that can be drawn from these figures is that chiral perturbation 
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theory for QCD to one-loop order is reasonably convergent throughout the PQ region defined 



by Xx < Xs 



X 



A, B,l,3. With our parameters, the least convergent quantity is fx- We 



also note that, generally speaking, the expansion is better behaved when one increases the 
masses along the lines at small angles. In other words, the expansion may be more reliable 
when valence masses are smaller than sea-quark masses, and vica-versa. This is good news 
for simulations, since pushing to small valence quark masses is relatively cheap. 

Of course, as noted in the previous section, the valence quark masses cannot be pushed 
too low at fixed sea masses. At some point the corrections 5*^ (in the vr-plane) and (in 
the K-plane) diverge. This is not apparent, however, from Fig. Bl. To observe the divergence 



we show in Fig. |T0| the region close to the x-axis, which does reveal the expected features. 
Though of theoretical interest, it seems that the breakdown of the chiral expansion due to 
enhanced chiral logarithms occurs only in a tiny region of parameter space and is therefore 
of little or no practical significance. 

Our final comment concerns the importance of including the non-analytic terms in fits 
to PQ data. The tree level contributions alone would have produced only straight lines in 
Figs. The curvature seen in the graphs is due to the logarithms originating in loop 
diagrams. An attempt to model data collected in the heavy quark mass region with only 
the tree level terms will clearly lead to a significant systematic errors in the extrapolations 
to QCD values. 



II 








0.2 0.4 0.6 0.8 



FIG. 10. Breakdown of the chiral expansion for small values of valence quark mass and a fixed 
sea quark mass. In the white region, the expansion breaks down as the NLO term in the expression 
for fx becomes greater than the LO term (5^ > !)• The following contours mark changes of .1 in 
6j^, so that the darkest region corresponds to 6^^ < .4 
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VI. DETERMINING L^ 

As noted above, the GL coefficients L4, L5, Lg and Lg can be obtained using PQ simu- 
lations with degenerate sea quarks. Simulations with non-degenerate sea quarks are useful 
but not essential. In this section we consider the other coefficient which enters into the 
NLO expressions for the physical meson masses, namely L7. This contributes only to flavor- 
diagonal mesons, i.e. the t] and ttq in QCD. Since the contribution is proportional to quark 
mass differences, the largest effect is on the t] mass. This is conveniently isolated using the 
violation of the mesonic Gell Mann-Okubo relation Wl\ 



3m^ log(m^) — log(m^) — Am\^ \og{m\, 
+ (6-1) 

Here we have given the result in the isospin-symmetric limit, i.e. with two degenerate quarks 
of mass m = {m^ + m^)/2, and a single strange quark.Q 

Thus one method for determining L7 is to calculate m^, mx and m.^ in unquenched 
simulations with non-degenerate sea quarks and ffi to the result in Eq. { p.l\) . In this section 
we show how one can obtain L7 using degenerate sea quarks by taking advantage of PQ 
simulations. 

A clue on how to proceed is provided by the fact that the rj propagator contains contrac- 
tions in which the quark propagators are disconnected ("hairpin" contractions). Close to 
the degenerate limit, these contractions are proportional to (m^ — m)^, due to cancelations 
between light and strange quark propagators. Comparing to Eq. (|6.1| ), and noting that 
ml — ml (X ms — m, it is plausible that the L7 contribution is related to the disconnected 
contraction, and that by studying this contraction in the PQ theory one can determine Lj. 
This is indeed what we find. 

First, we recall the general form of the tree level propagators. In section |T| we showed 
that the propagators for flavor off-diagonal mesons, G^^, have an ordinary single pole at 
the meson mass [Eq. (|3.1|)], while G^, the "neutral" or flavor-diagonal propagators, have a 
double pole at the same mass if a = 6 or ma = mt, [Eq. (|3.3| )]. As we show below, these 
general forms prevail also at NLO: 

Gab{p) \mA=n.B = „2 , + ^ou-pole , (6.2) 



Z T) 

G'ab(p) L^=ms = (p2 + M2^^2 + single-pole + non-pole , (6.3) 



Z T) 

— — ^ ^ + single-pole -|- non-pole . (6.4) 



Note that the one-loop expressions for the dependence of txik and m,r on quark masses can be 
obtained from the general results above by setting A^i = = -^3 = !> = "t-I) and choosing the 
valence masses to be equal to the appropriate sea quark masses. 
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Here Maa is a shorthand for MAsi'mA = tt^b), and its NLO expression is given in Eq. ( [4.1D . 
Note that these equations are vahd only if A and B are partially quenched rather than 
unquenched (i.e. ttia 7^ mi), and apply only close to the poles at p"^ = M\^. In particular, 
the neutral propagators have additional poles, at the masses of the neutral sea-quark mesons. 
These are not of interest here. 

The quantity we propose to use to determine L7 is V, the (suitably normalized) coefficient 
of the double pole in the neutral propagators. An equivalent definition in terms of lattice 
observables is 



J d^x{TrAA{t,x)7!-BB{0)) 



— > vt , , 

. (6.5) 

mA=mB t ^ 00 2MaA 



J d^x{TTABit,x)TTBAiO)) 

Here the fiavor indices are chosen to select the "disconnected" (numerator) and "connected" 
(denominator) contractions contributing to the ttaa propagator .0 At large times, the denom- 
inator is dominated by the single-pole contribution of the pseudo-Goldstone boson of mass 
Maa- Only double-pole contributions to the numerator lead to the ratio growing linearly 
with t, and V measures their size. We note also that our ratio is the standard one used in 
studies of artifacts in the quenched theory. 

We claim that P is a "physical" quantity in the PQ theory, on a par with the valence 
meson masses such as Maa- In support of this claim, we note that V is determined from 
the long distance properties of correlators, and is independent of the choice of interpolating 
fields (because it is determined by a ratio). In particular, we claim that T> can be calculated 
using the effective low-energy theory with the same level of reliability as the meson masses.[] 
We stress that the coefficients of the single poles, unlike the double pole, do depend on 
the choice of interpolating fields and are not quantities which can be predicted using chiral 
perturbation theory. 

In PQChPT, we obtain V by calculating the propagators G^) and at NLO and using 
Eqs. ( |6.2| )-( |6l4D . The tree level result for V can be read off from Eqs. ( p.l|) and (|3] 



-p _ 1 (Xi - Xa) (X2 - Xa) (X3 - Xa) ,ggx 
N iX7r-XA){Xv-XA) 

Note that the residue of the double-pole vanishes whenever the valence quark mass equals 
any of the sea quark masses. This must be the case since one is then considering a correlator 
which could be constructed entirely from sea quarks, and thus is physical, and cannot contain 
double poles. 

The dependence of V on Ly begins at one-loop order. We have calculated V to this order 
only for the case of degenerate sea quarks (mi = m2 = m^). Details are given in app. ^. 
The result is 



^ This choice simplifies the numerical calculation, but is not necessary. It follows from Eq. ( |6.4D 
that the numerator could be replaced by {7rAA{t,p = 0)ttaa{0)), i.e. the temporal Fourier transform 
of G^y^. This changes the t— independent part of the ratio but leaves V unchanged. 

^Proving this claim rigorously seems difficult, given the lack of a Hamiltonian formulation of the 
PQ theory. Our point, however, is that the calculation of P is as well controlled as that of meson 
masses. 
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V = (ML - Mi A -^\L' — ^ + — (Mia - Mi, 

+ ^^^2^2 (^[^^ ~ XAf\ogxsA + x\^og{xA/xsA) + xl^og{xs / Xsa)^ , (6.7) 

where we use S to denote the sea-quark, and xsa = (X5 + Xa)/2. The first term is the same 
as the tree-level result with one-loop corrected meson mass-squareds replacing quark masses. 
The second term is the analytic term containing the L'j dependence. As advertised, L'j and 
TrtQ appear in the appropriate combination to be combined into the standard L7 [Eq. (|2.14| )]. 



The logarithmic terms, from wavefunction renormalization and from loop diagrams, combine 
into a fairly simple form. One check on the result is that the anomalous dimension of L5 is 
such that it cancels the dependence on the choice of the scale in the logarithm (for = 3, 
where the anomalous dimension is known). Another is that it vanishes when mi = tjia- 

Thus, from the coefficient of the double pole, one can extract the combination 2NL'! + L^. 
Combined with the results of the previous sections this allows a determination of L7. 



VII. CONCLUSIONS 

Partially quenched theories can play an important role in determining physical parame- 
ters of QCD. Our results show how one can use them to simplify the extrapolation to QCD 
and the determination of the GL parameters L4_8. In particular, it is sufficient, though not 
necessary, to use degenerate sea quarks. One must, however, use three sea quarks. 

Our approach relies on chiral perturbation theory at next-to-leading order. An important 
issue is how light the sea quarks need to be for our formulae to be sufficiently accurate. A 
conservative approach is to work down to masses where the NLO corrections themselves are 



10% of the leading order result, so that the missing NNLO terms are very small [0. This 
requires working down to mgea ~ m^/S. Another approach is to look at our figures and 
see what mass range is required to observe the predicted curvature. To do so appears to 
require working down to at least mgea ~ In the end, this issue can be resolved using 

simulations themselves, including partially quenched simulations, to check the reliability of 
the NLO predictions. 

This exercise will also shed light on the question of whether the physical strange quark 
mass is light enough that NLO chiral perturbation theory is applicable for QCD. Even if 
the strange quark turns out to be too heavy, the results presented here still apply to the 
light-quark sector, if we set N = 2. 

It is often observed that linear fits are adequate for the mass dependence of physical 
quantities in the range mgea = ^s/'^ — ttLs. Our results make clear, however, that extending 
such linear fits to the chiral limit, i.e. leaving out the curvature due to chiral logarithms, 
can lead to a significant error in the determination of physical parameters. This is most 
clearly illustrated by Fig. |. 

All the formulae we present are for infinite volume, and thus apply for box sizes such that 
Mt^L ^ 1. As pion masses decrease, this requires working in boxes of increasing size. This 
should not, however, be necessary. Chiral perturbation theory remains valid for finite Mt^L, 
as long as Aqcd-^ ^ 1, and can be used to calculate the volume dependence of physical 
quantities. It would be interesting to do this for the quantities we consider here. 
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Another interesting extension of our calculations is to include the effect of discretiza- 
tion errors within the chiral Lagrangian itself, i.e. to use the appropriate Lagrangian for 
the lattice theory at non-zero lattice spacing. This would aid in the extrapolation to the 
continuum limit. 

Finally, we comment on the theoretical status of our methods. We use a Lagrangian 
containing the v( (or more precisely the "super-?]'", $o); and our calculations rely on as- 
sumptions about the size of its couplings. One can show, however, that there is a limit in 
which our calculation is equivalent to that with the r( integrated out non-perturbatively . 
In that limit (essentially mg — > oo) the assumptions that we have made are valid. One might 
also be concerned about the theoretical foundation of the whole calculation: Is it justified 
to use chiral perturbation theory for the unphysical PQ theory? We have also made some 
progress on this issue. One can show that, in certain cases, derivatives of partially quenched 
quantities with respect to valence quark masses can be exactly related to derivatives of un- 
quenched quantities with respect to sea quark masses [Q. Our NLO results are consistent 
with these exact relations. 
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APPENDIX A: NEUTRAL MESON PROPAGATOR AT TREE LEVEL 

Here we calculate the neutral meson propagator. From £lo we find 

Gjv' = Go 1 + V (Al) 

ipt)^, =(P' + Xa)5afeea (A2) 

V..= '^^e.e,. (A3) 

The full propagator is thus 

= (Go ^ + vy' = {l + GoVy'Go . (A4) 
Because V is an outer product, the combination GqV is proportional to a projection operator 

GqV 
tr(GoV^) ' 



A^-^^, A' = A. (A5) 



Thus for any function /, 

fiA) - m = A[f{l) - /(O)] , (A6) 

and so 
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1 



A 



1 + tT{GoV) 
GoV 



Inserting this in Eq. 



we find 

Gn = Gq 



tr(GoV) ■ 



GqVGo 
1 + tr(GoV^) ■ 



(A7) 



(A8) 



which reproduces the result of Ref . . 

The analytic structure of the propagator is not clear from this result. In particular, 
its diagonal elements [GAr]aa appear to contain double poles (from the two factors of Gq in 
the second term). We know, however, that if we restrict ourselves to the physical sea-quark 
sector then Gn cannot contain double poles. Thus there are cancelations hidden in Eq. ( |A8| ) 
which we want to make explicit. 

To do so we need to introduce the restriction of the various matrices to the sea-sea 
sector. We denote these restrictions by overbars. We first observe that the previous steps 
go through identically for the restricted matrices 



{Gn 



-1 



Gn 



Gn V Gn 



1 + tr(GoF) ' 



(A9) 



Comparing this with the restriction of Eq. (|A8| 



Gm — Gn 



GqVGq 

1 + tr(Go^) ' 



and using the result 



tr(Gol^) = tr(G'oV^) , 
(which follows because the valence and ghost contributions cancel) we find 

Gn = [Gn^ 



-1 



(AlO) 



(All) 



(A12) 



In other words, restriction to the sea-quark sector commutes with inversion. This result is 
non-trivial because is not block diagonal — it encapsulates the lack of feedback from the 
valence to the sea-quark sector. 

Returning to the simplification of the propagator, we note that 



det 



N 



G 



/det 



Gq 



det 
det 



Go Gj^ 



l + tr{GoV)A 



exptrln l + ti(GoV)A 
expln [l +tr(GoF) 
1 + tr(GoF) . 



(A13) 
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Thus the factor multiplying the double pole in Eq. ([A8|) can be rewritten as 



[l + tr(G'o1/)]"^ = det 
The determinant in the numerator is simple 



/det 



G 



(AM) 



det 



\N3 



(A15) 



Our final task is to evaluate the determinant in the denominator. 

To do this, we note that = G]^^ is block diagonal. The exact SU{Ni) x SU{N2) x 
SUlN^) fiavor symmetry implies that there are, for each sea-quark type i, Ni — 1 fiavor 
non-singlet neutral pions which are eigenvectors of Gjy • Since V projects onto fiavor singlet 
states, the corresponding eigenvalues are those of Gq , namely (p^ + Xi)- The non-trivial 
part of Gj^ , to which V does contribute, is thus a 3 x 3 block. In QCD, with Ni = 1, this is 
the entire matrix, and describes the ttq, rj and 1]'. For convenience, we use these names for 
general iVj as well. We denote the restriction of matrices to this subspace with double bars. 
A straightforward exercise shows that 



SR^DRS . 



S = diag(l, 1, a/1 + a) , 

D = diag(/ + X7r,P^ + + Xv' 



(A16) 
(A17) 
(A18) 



with R a rotation matrix. S rescales the singlet field so that it has a canonical kinetic term. 
The meson mass-squareds are given, up to corrections of size x^ /Tn^, by 



Xr,' = ("^0 + X)/(l + «) 



(A19) 
(A20) 
(A21) 



where 



X 



AT X! ^iXi 



N 



X~ 



i=l,3 



N 



4 = 1,3 



(A22) 



are averages over the sea sector. 
Using these results we find that 



det 



N 



G 



(p' + Xi)^^-^(p' + X2)^^^-^(p^ + X3)^^-'(p' + X.)(p^ + X.)(p^ + X.')(l+«)- (A23) 



Inserting this and Eq. ( |A15| ) into Eqs. ( |A14| ) and ( [A8| ) gives the result quoted in the main 
text, Eq. ( p.3|) . It is straightforward to generalize this result to an arbitrary number of 
different sea-quark masses. 
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APPENDIX B: ONE-LOOP CALCULATION OF V 



To extract V using Eqs. ( |6.4|) and ( |6.2|) we need the NLO results for the charged and 
neutral propagators. We do the calculation only for N degenerate sea quarks, which we 
denote using the label S rather than 1, 2, 3. 

The NLO calculation of the charged propagator was described in Sec. fV], and the ex- 
pression for M\j^ can be obtained using Eq. ( [4.1|) . Here we also need the wavefunction 
renormalization factor, for which we find 



T ^ (T AT , r A ^ 1 XS + XA. (XS + XA 

1 - —[L^Nxs + L^Xa) - -ITTTT^ ^ log 



(Bl) 



For the neutral propagator we need to generalize the calculation of app. |^ to NLO. The 
inverse propagator becomes 

G~N = Go ' + 1^ + S , (B2) 

where S contains the NLO tree and one-loop contributions. The structure of S is such that 
the method used to calculate in Sec. |l|does not apply, and we simply invert the matrix 
by brute force. 

For our purposes, it is sufficient to consider the restriction of Ga? to the three-dimensional 
basis 

1 ^ 

This is because, first, we want only the tiaa propagator and so do not need to introduce an 
additional valence quark; and, second, because we use degenerate sea quarks so that there 
is no mixing of t^aa with fiavor non-singlet neutral sea-quark mesons. The generalization to 
non-degenerate sea quarks, which involves a larger basis of neutral states, is straightforward 
in principle, but tedious in practice, and we have not carried it out. 

The contributions to S fall into two classes: those that are common to the charged mesons 
Sc, and those that are special to the neutral mesons Stv- The former can be obtained from 
the results of Sec. |^ and we find 

Go ^ + Sc = diag(t;, w, -v) , (B4) 
v={p' + MI^)/Za, (B5) 
w={/ + Mls)/Zs. (B6) 

Here M'^g is the squared mass of fiavor non-singlet mesons composed of sea quarks evaluated 
at NLO, and Zs the corresponding wavefunction renormalization. Their expressions can be 
obtained from those for M^^ and Za by the substitution Xa^ Xs- 

We now turn to S^r. Because of the graded symmetry, it shares a particular matrix 
structure with V , and it is convenient to make this explicit by the following definitions: 




1/ + Ejv = ^/Ny Nz -VNy . (B7) 
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The contributions to this matrix from V are 

xv = Vv = zv = xo = {ml + ap'^)/N . (B8) 
Thus we write x = Xq + 6x, y = Xq + 6y and z = xq + 6z. 




NLO LO 



FIG. 11. Diagrams contributing to Sjv. The letters next to the hnes denote the flavor indices 
of the propagating mesons. "LO" and "NLO" describe the order of the vertex that makes the 
diagram contribute at NLO. 

The diagrams contributing to Sjy are of the general form shown in Fig. |TT]. The tree 
level contributions come from the two- meson vertices included in Cnlo,2 [Eq. ( |2.7|) ]. The V2 
term acts as a subleading correction to mg. Since at LO [Eq. ( |6.6| )], V is independent of ml, 
it follows that at NLO V can at most depend on its leading value. Thus V2 contributes only 
at NNLO. The other tree level diagram comes from the L7 term which gives the following 
contributions: 

SW^ = l-2Lrx\, = ^2L7XAX5 , S^'h = ^2L,xl ■ (B9) 



FIG. 12. Each line in the diagram represents a single quark field. This is the only loop diagram 
of the type that appears in Fig. 11 in which the quark lines are disconnected. 



Figure |T2| shows the quark line structure of the only one-loop graph involving the lowest 
order vertices which corresponds to disconnected quark lines, and so does not contribute to 
Ec. The contribution to Sat is 



6^'^x 
5(2);, 



1 



(p^ - 2xa)xa logXA 



167r2/^3 



olP -Xa-Xs) 7; log 



Xa + Xs 



\{P^ - '^Xs)Xs^ogxs ■ 



167r2/2 3 



(BIO) 

(Bll) 
(B12) 
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Collecting these contributions we end up with 



V + X V Ny 
G-i = I w + Nz 

—X —y/Ny 




(B13) 



The relevant part of the inverse is 



Gaa = - 



1 1 xw + N{xz — y"^) 



V 



w + Nz 



(B14) 



The first term is the one-loop corrected single pole, while the second contains the expected 
double pole. Inserting this result into the definitions Eqs. ([61^ ) and ( |6.4|) , and using Eqs. ( pi]) 
and (|B5D, we can read off the required double-pole coefficient 



Expanding in powers of x we find 



xw + N{xz — y^) 
w + Nz 



(B15) 



(B16) 



Note that {xyZy — yy) = 0, so that a possible contribution proportional to ml cancels. 
Substituting and rearranging, we find the answer Eq. (|6l7|). 
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